Genomic insights into anthropozoonotic tuberculosis in captive sun bears (Helarctos malayanus) and an Asiatic black bear (Ursus thibetanus) in Cambodia.

Contact between humans and wildlife presents a risk for both zoonotic and anthropozoonotic disease transmission. In this study we report the detection of human strains of Mycobacterium tuberculosis in sun bears and an Asiatic black bear in a wildlife rescue centre in Cambodia, confirming for the first time the susceptibility of these bear species to tuberculosis when in close contact with humans. After genotyping revealed two different strains of M. tuberculosis from cases occurring between 2009 and 2019, 100 isolates from 30 sun bear cases, a single Asiatic black bear case, and a human case were subjected to whole genome sequencing. We combined single nucleotide polymorphism analysis and exploration of mixed base calls with epidemiological data to indicate the evolution of each outbreak. Our results confirmed two concurrent yet separate tuberculosis outbreaks and established a likely transmission route in one outbreak where the human case acted as an intermediatory between bear cases. In both outbreaks, we observed high rates of transmission and progression to active disease, suggesting that sun bears are highly susceptible to tuberculosis if exposed under these conditions. Overall, our findings highlight the risk of bi-directional transmission of tuberculosis between humans and captive bears in high human tuberculosis burden regions, with implied considerations for veterinary and public health. We also demonstrate the use of standard genomic approaches to better understand disease outbreaks in captive wildlife settings and to inform control and prevention measures.

Taking a color photo: A homozygous 25-bp deletion in Bace2 may cause brown-and-white coat color in giant pandas.

Brown-and-white giant pandas (hereafter brown pandas) are distinct coat color mutants found exclusively in the Qinling Mountains, Shaanxi, China. However, its genetic mechanism has remained unclear since their discovery in 1985. Here, we identified the genetic basis for this coat color variation using a combination of field ecological data, population genomic data, and a CRISPR-Cas9 knockout mouse model. We de novo assembled a long-read-based giant panda genome and resequenced the genomes of 35 giant pandas, including two brown pandas and two family trios associated with a brown panda. We identified a homozygous 25-bp deletion in the first exon of Bace2, a gene encoding amyloid precursor protein cleaving enzyme, as the most likely genetic basis for brown-and-white coat color. This deletion was further validated using PCR and Sanger sequencing of another 192 black giant pandas and CRISPR-Cas9 edited knockout mice. Our investigation revealed that this mutation reduced the number and size of melanosomes of the hairs in knockout mice and possibly in the brown panda, further leading to the hypopigmentation. These findings provide unique insights into the genetic basis of coat color variation in wild animals.

Characterization of antibiotic resistance genes and mobile genetic elements in Escherichia coli isolated from captive black bears.

The objective of this study was to analyze the antimicrobial resistance (AMR) characteristics produced by antibiotic resistance genes (ARGs), mobile genetic elements (MGEs) and gene cassettes in Escherichia coli isolated from the feces of captive black bears. Antimicrobial susceptibility testing was performed by using the disk diffusion method, and both MGEs and integron gene cassettes were detected by polymerase chain reaction. Our results showed that 43.7% (62/142) of the isolates were multidrug resistant strains and 97.9% (139/142) of the isolates were resistant to at least one antibiotic. The highest AMR phenotype was observed for tetracycline (79.6%, 113/142), followed by ampicillin (50.0%, 71/142), trimethoprim-sulfamethoxazole (43.7%, 62/142) and cefotaxime (35.9%, 51/142). However, all isolates were susceptible to tobramycin. tetA had the highest occurrence in 6 ARGs in 142 E. coli isolates (76.8%, 109/142). Ten mobile genetic elements were observed and IS26 was dominant (88.0%, 125/142). ISECP1 was positively associated with five ß-lactam antibiotics. ISCR3/14, IS1133 and intI3 were not detected. Seventy-five E. coli isolates (65 intI1-positive isolates, 2 intI2-positive isolates and 8 intI1 + intI2-positive isolates) carried integrons. Five gene cassettes (dfrA1, aadA2, dfrA17-aadA5, aadA2-dfrA12 and dfrA1-aadA1) were identified in the intI1-positive isolates and 2 gene cassettes (dfrA1-catB2-sat2-aadA1 and dfrA1-catB2-sat1-aadA1) were observed in the intI2-positive isolates. Monitoring of ARGs, MGEs and gene cassettes is important to understand the prevalence of AMR, which may help to introduce measures to prevent and control of AMR in E. coli for captive black bears.

Evidence of population genetic structure in Ecuadorian Andean bears.

Wildlife conservation in Andean countries is a global priority because of the high levels of biodiversity and endemism. Historically, these countries have had limited resources to monitor wildlife (e.g., through genetic tools) and establish conservation programs. Focusing on the study and emblematic use of a few charismatic species has been a strategic approach to direct efforts for conservation and development planning. Consequently, the Andean bear is a flagship and umbrella species for highly biodiverse Andean countries like Ecuador. The few studies exploring the population genetics of this species have concluded that it has low genetic diversity and few units for conservation as populations appear to be well connected. However, these results might be attributed to ascertainment bias as studies have been performed with heterologous molecular markers. Here, using both mtDNA sequences and species-specific microsatellite markers, we show that Andean bears in Ecuador have population structure. Additionally, we found through the study of three Ecuadorian populations that the species might have a higher genetic diversity than we previously thought. These results could support the revision of research priorities, conservation, and planning strategies to improve connectivity for this species which occurs in crucial biodiversity hotspots.

Single-cell mRNA sequencing of giant panda (Ailuropoda melanoleuca) seminoma reveals the cellular and molecular characteristics of tumour cells.

Testicular tumours are zoonoses that can occur in not only human, but other animals, include giant pandas. A middle-aged male giant panda named Fufu was diagnosed with a testicular tumour and underwent surgery to remove the entire left testis. The testis was mainly composed of three substantive parts: normal tissue on the outside, tumour tissue in the middle, and necrosis in the centre. HE stains revealed that the tumour was a seminoma. Single-cell mRNA sequence was applied to characterise cellular states and molecular circuitries of giant panda testicular seminoma. Only germ cell markers expressed in nearly all tumour cells, and the tumour cells appeared to be the same subtype of seminoma cells. We identified four clusters with unique genes expression. They were early apoptosis cells (EAC), inactive cells (IC), active cells subcluster 1 (AC-1) and active cells subcluster 2 (AC-2). We utilised monocle tools and found that IC cells was in the initiation stage, and EAC was one type of terminal stage, suggesting that tumour cells may undergo apoptosis in the future. AC-2 was another type of terminal stage, representing a group of progressive cells. Our study represents the first report to utilise scRNA-seq to characterise the cellular states and molecular circuitries of a giant panda testicular tumour. This investigation proposes CD117 and CD30 as dependable markers for future pathologic diagnosis. Our findings also suggest that CTSV and other genes with unique expression patterns in active and progressive giant panda seminoma cells may act as early prognostic biomarkers.

Evaluating hybrid speciation and swamping in wild carnivores with a decision-tree approach.

Hybridization is an important evolutionary force with a principal role in the origin of new species, known as hybrid speciation. However, ongoing hybridization can create hybrid swamping, in which parental genomes are completely lost. This can become a biodiversity threat if it involves species that have adapted to certain environmental conditions and occur nowhere else. Because conservation scientists commonly have a negative attitude toward hybrids, it is important to improve understanding of the influence of interspecific gene flow on the persistence of species. We reviewed the literature on species hybridization to build a list of all known cases in the order Carnivora. To examine the relative impact, we also noted level of introgression, whether fertile offspring were produced, and whether there was mention of negative or positive evolutionary effects (hybrid speciation and swamping). To evaluate the conservation implications of hybrids, we developed a decision-making tree with which to determine which actions should be taken to manage hybrid species. We found 53 hybrids involving 68 unique taxa, which is roughly 23% of all carnivore species. They mainly involved monophyletic (83%) and sympatric species (75%). For 2 species, the outcome of the assessment was to eliminate or restrict the hybrids: Ethiopian wolf (Canis simensis) and Scottish wildcat (Felis silvestris silvestris). Both species hybridize with their domestic conspecifics. For all other cases, we suggest hybrids be protected in the same manner as native species. We found no evidence of genomic extinction in Carnivora. To the contrary, some species appear to be of hybrid origin, such as the Asiatic black bear (Ursus thibetanus) and African golden wolf (Canis lupaster). Other positive outcomes of hybridization are novel genetic diversity, adaptation to extreme environments, and increased reproductive fitness. These outcomes are particularly valuable for counterbalancing genetic drift and enabling adaptive introgression in a human-dominated world.

La especiación por hibridación es una fuerza evolutiva importante con un papel principal en el origen de una nueva especie. Sin embargo, la hibridación continua puede generar un estancamiento híbrido en el que se pierden por completo los genomas parentales. Esto puede convertirse en una amenaza para la biodiversidad si involucra a una especie que se ha adaptado a ciertas condiciones ambientales y sólo se encuentra en un lugar. Ya que los científicos de la conservación suelen tener una actitud negativa hacia los híbridos, es importante incrementar el entendimiento de la influencia que tiene el flujo interespecífico sobre la persistencia de las especies. Revisamos la literatura sobre la hibridación de especies para generar una lista de todos los casos conocidos en el orden Carnívora. También observamos el nivel de introgresión, si se produjo descendencia fértil y si hubo mención de los efectos evolutivos positivos o negativos (especiación híbrida y estancamiento) para analizar el impacto relativo. Desarrollamos un árbol de decisión con el cual determinar cuáles acciones deberían tomarse en el manejo de las especies híbridas para evaluar las implicaciones que tienen los híbridos para la conservación. Encontramos 53 híbridos de 68 taxones únicos, lo que representa aproximadamente el 23% de todos los carnívoros. Estos híbridos incluyen principalmente a especies monofiléticas (83%) y simpátricas (75%). Para dos especies, los resultados del análisis fueron la eliminación o restricción de los híbridos: el lobo etíope (Canis simensis) y el lince escocés (Felis silvestris silvestris). Ambas especies hibridan con sus coespecíficos domésticos. Para todos los demás casos sugerimos que se proteja a los híbridos de la misma manera que a las especies nativas. No encontramos evidencias de una extinción genómica en el orden Carnívora. Al contrario, algunas especies parecen tener un origen híbrido, como el oso negro asiático (Ursus thibetanus) y el lobo dorado africano (Canis lupaster). Otros resultados positivos de la hibridación son la diversidad genética novedosa, la adaptación a ambientes extremos y el incremento en la adaptabilidad reproductiva. Estos resultados son de valor particular para contrarrestar la deriva génica y permitir la introgresión adaptativa en un mundo dominado por humanos. Evaluación de la especiación y estancamiento en carnívoros silvestres con una estrategia de árbol de decisión.

Conservation Genomics and Metagenomics of Giant and Red Pandas in the Wild.

Giant pandas and red pandas are endangered species with similar specialized bamboo diet and partial sympatric distribution in China. Over the last two decades, the rapid development of genomics and metagenomics research on these species has enriched our knowledge of their biology, ecology, physiology, genetics, and evolution, which is crucial and useful for their conservation. We describe the evolutionary history, endangerment processes, genetic diversity, and population structure of wild giant pandas and two species of red pandas (Chinese and Himalayan red pandas). In addition, we explore how genomics and metagenomics studies have provided insight into the convergent adaptation of pandas to the specialized bamboo diet. Finally, we discuss how these findings are applied to effective conservation management of giant and red pandas in the wild and in captivity to promote the long-term persistence of these species.

Individual genotypes from environmental DNA: Fingerprinting snow tracks of three large carnivore species.

Continued advancements in environmental DNA (eDNA) research have made it possible to access intraspecific variation from eDNA samples, opening new opportunities to expand non-invasive genetic studies of wildlife populations. However, the use of eDNA samples for individual genotyping, as typically performed in non-invasive genetics, still remains elusive. We present successful individual genotyping of eDNA obtained from snow tracks of three large carnivores: brown bear (Ursus arctos), European lynx (Lynx lynx) and wolf (Canis lupus). DNA was extracted using a protocol for isolating water eDNA and genotyped using amplicon sequencing of short tandem repeats (STR), and for brown bear a sex marker, on a high-throughput sequencing platform. Individual genotypes were obtained for all species, but genotyping performance differed among samples and species. The proportion of samples genotyped to individuals was higher for brown bear (5/7) and wolf (7/10) than for lynx (4/9), and locus genotyping success was greater for brown bear (0.88). The sex marker was typed in six out of seven brown bear samples. Results for three species show that reliable individual genotyping, including sex identification, is now possible from eDNA in snow tracks, underlining its vast potential to complement the non-invasive genetic methods used for wildlife. To fully leverage the application of snow track eDNA, improved understanding of the ideal species- and site-specific sampling conditions, as well as laboratory methods promoting genotyping success, is needed. This will also inform efforts to retrieve and type nuclear DNA from other eDNA samples, thereby advancing eDNA-based individual and population-level studies.

Single-nucleus transcriptome inventory of giant panda reveals cellular basis for fitness optimization under low metabolism.

BACKGROUND: Energy homeostasis is essential for the adaptation of animals to their environment and some wild animals keep low metabolism adaptive to their low-nutrient dietary supply. Giant panda is such a typical low-metabolic mammal exhibiting species specialization of extremely low daily energy expenditure. It has low levels of basal metabolic rate, thyroid hormone, and physical activities, whereas the cellular bases of its low metabolic adaptation remain rarely explored. RESULTS: In this study, we generate a single-nucleus transcriptome atlas of 21 organs/tissues from a female giant panda. We focused on the central metabolic organ (liver) and dissected cellular metabolic status by cross-species comparison. Adaptive expression mode (i.e., AMPK related) was prominently displayed in the hepatocyte of giant panda. In the highest energy-consuming organ, the heart, we found a possibly optimized utilization of fatty acid. Detailed cell subtype annotation of endothelial cells showed the uterine-specific deficiency of blood vascular subclasses, indicating a potential adaptation for a low reproductive energy expenditure. CONCLUSIONS: Our findings shed light on the possible cellular basis and transcriptomic regulatory clues for the low metabolism in giant pandas and helped to understand physiological adaptation response to nutrient stress.

The bear necessities: A sensitive qPCR assay for bear DNA detection from bile and derived products to complement wildlife forensic enforcement.

Demand for bear bile, a prized component of traditional Asian medicines, threaten Asiatic and sun bear population sustainability. While laws exist to prevent poaching and trafficking of bear parts and derivatives, smuggling persists with demand extending to surrogate species, including American black bears (Ursus americanus). Mitochondrial DNA (mtDNA) sequencing can identify products putatively containing biological bear material but can be undermined by PCR inhibitors in bile and a lack of sensitivity at trace levels. Quantitative PCR (qPCR) assays can be used to distinguish between closely related target species, while concomitantly evaluating inhibition and false negative results in low quality/quantity DNA applications. Herein, we develop a multiplexed qPCR assay to detect and differentiate among bear species, including highly diluted bile samples mixed within liquors as common dilutants. The assay detects as little as 10 locus copies/reaction of bear DNA with 95% confidence, distinguishing among sun, Asiatic and American black bears. Demonstrating the sensitivity and applicability of this assay in context of current bile mixture recipes, dilutions of 1:5,000 bile with ethanol, red wine, and spirits, all yielded clear quantifiable detections, where our data suggests as little as 1 drop of bile per 750 mL bottle of alcohol would still exceed the limits of detection (e.g., 1:15000 dilution or <0.05 mL bile per 750 mL bottle). Overall, this study provides a rapid, sensitive, and specific test to identify and distinguish among bear species commonly used for bile production to aid wildlife enforcement applications.

Circadian gene transcription plays a role in cellular metabolism in hibernating brown bears, Ursus arctos.

Hibernation is a highly seasonal physiological adaptation that allows brown bears (Ursus arctos) to survive extended periods of low food availability. Similarly, daily or circadian rhythms conserve energy by coordinating body processes to optimally match the environmental light/dark cycle. Brown bears express circadian rhythms in vivo and their cells do in vitro throughout the year, suggesting that these rhythms may play important roles during periods of negative energy balance. Here, we use time-series analysis of RNA sequencing data and timed measurements of ATP production in adipose-derived fibroblasts from active and hibernation seasons under two temperature conditions to confirm that rhythmicity was present. Culture temperature matching that of hibernation body temperature (34 °C) resulted in a delay of daily peak ATP production in comparison with active season body temperatures (37 °C). The timing of peaks of mitochondrial gene transcription was altered as were the amplitudes of transcripts coding for enzymes of the electron transport chain. Additionally, we observed changes in mean expression and timing of key metabolic genes such as SIRT1 and AMPK which are linked to the circadian system and energy balance. The amplitudes of several circadian gene transcripts were also reduced. These results reveal a link between energy conservation and a functioning circadian system in hibernation.

Range-wide evolutionary relationships and historical demography of brown bears (Ursus arctos) revealed by whole-genome sequencing of isolated central Asian populations.

Phylogeographic studies uncover hidden pathways of divergence and inform conservation. Brown bears (Ursus arctos) have one of the broadest distributions of all land mammals, ranging from Eurasia to North America, and are an important model for evolutionary studies. Although several whole genomes were available for individuals from North America, Europe and Asia, limited whole-genome data were available from Central Asia, including the highly imperilled brown bears in the Gobi Desert. To fill this knowledge gap, we sequenced whole genomes from nine Asian brown bears from the Gobi Desert of Mongolia, Northern Mongolia and the Himalayas of Pakistan. We combined these data with published brown bear sequences from Europe, Asia and North America, as well as other bear species. Our goals were to determine the evolutionary relationships among brown bear populations worldwide, their genetic diversity and their historical demography. Our analyses revealed five major lineages of brown bears based on a filtered set of 684,081 single nucleotide polymorphisms. We found distinct evolutionary lineages of brown bears in the Gobi, Himalayas, northern Mongolia, Europe and North America. The lowest level of genetic diversity and the highest level of inbreeding were found in Pakistan, the Gobi Desert and Central Italy. Furthermore, the effective population size (Ne ) for all brown bears decreased over the last 70,000 years. Our results confirm the genetic distinctiveness and ancient lineage of brown bear subspecies in the Gobi Desert of Mongolia and the Himalayas of Pakistan and highlight their importance for conservation.

Corridor-based approach with spatial cross-validation reveals scale-dependent effects of geographic distance, human footprint and canopy cover on grizzly bear genetic connectivity.

Understanding how human infrastructure and other landscape attributes affect genetic differentiation in animals is an important step for identifying and maintaining dispersal corridors for these species. We built upon recent advances in the field of landscape genetics by using an individual-based and multiscale approach to predict landscape-level genetic connectivity for grizzly bears (Ursus arctos) across ~100,000 km2 in Canada's southern Rocky Mountains. We used a genetic dataset with 1156 unique individuals genotyped at nine microsatellite loci to identify landscape characteristics that influence grizzly bear gene flow at multiple spatial scales and map predicted genetic connectivity through a matrix of rugged terrain, large protected areas, highways and a growing human footprint. Our corridor-based modelling approach used a machine learning algorithm that objectively parameterized landscape resistance, incorporated spatial cross validation and variable selection and explicitly accounted for isolation by distance. This approach avoided overfitting, discarded variables that did not improve model performance across withheld test datasets and spatial predictive capacity compared to random cross-validation. We found that across all spatial scales, geographic distance explained more variation in genetic differentiation in grizzly bears than landscape variables. Human footprint inhibited connectivity across all spatial scales, while open canopies inhibited connectivity at the broadest spatial scale. Our results highlight the negative effect of human footprint on genetic connectivity, provide strong evidence for using spatial cross-validation in landscape genetics analyses and show that multiscale analyses provide additional information on how landscape variables affect genetic differentiation.

Characteristics of mesenchymal stem cells and their exosomes derived from giant panda (Ailuropoda melanoleuca) endometrium.

Conservation of genetic resources is an important way to protect endangered species. At present, mesenchymal stem cells (MSCs) have been isolated from the bone marrow and umbilical cords of giant pandas. However, the types and quantities of preserved cell resources were rare and limited, and none of MSCs was derived from female reproductive organs. Here, we first isolated MSCs from the endometrium of giant panda. These cells showed fibroblast morphology and expressed Sox2, Klf4, Thy1, CD73, CD105, CD44, CD49f, and CD105. Endometrium mesenchymal stem cells (eMSCs) of giant panda could induce differentiation into three germ layers in vitro. RNA-seq analysis showed that 833 genes were upregulated and 716 genes were downregulated in eMSCs compared with skin fibroblast cells. The results of GO and the KEGG analysis of differentially expressed genes (DEGs) were mainly focused on transporter activity, signal transducer activity, pathways regulating pluripotency of stem cells, MAPK signaling pathway, and PI3K-Akt signaling pathway. The genes PLCG2, FRK, JAK3, LYN, PIK3CB, JAK2, CBLB, and MET were identified as hub genes by PPI network analysis. In addition, the exosomes of eMSCs were also isolated and identified. The average diameter of exosomes was 74.26 ± 13.75 nm and highly expressed TSG101 and CD9 but did not express CALNEXIN. A total of 277 miRNAs were detected in the exosomes; the highest expression of miRNA was the has-miR-21-5p. A total of 14461 target genes of the whole miRNAs were predicted and proceeded with functional analysis. In conclusion, we successfully isolated and characterized the giant panda eMSCs and their exosomes, and analyzed their functions through bioinformatics techniques. It not only enriched the conservation types of giant panda cell resources and promoted the protection of genetic diversity, but also laid a foundation for the application of eMSCs and exosomes in the disease treatment of giant pandas.

Feeding during hibernation shifts gene expression toward active season levels in brown bears (Ursus arctos).

Hibernation in bears involves a suite of metabolical and physiological changes, including the onset of insulin resistance, that are driven in part by sweeping changes in gene expression in multiple tissues. Feeding bears glucose during hibernation partially restores active season physiological phenotypes, including partial resensitization to insulin, but the molecular mechanisms underlying this transition remain poorly understood. Here, we analyze tissue-level gene expression in adipose, liver, and muscle to identify genes that respond to midhibernation glucose feeding and thus potentially drive postfeeding metabolical and physiological shifts. We show that midhibernation feeding stimulates differential expression in all analyzed tissues of hibernating bears and that a subset of these genes responds specifically by shifting expression toward levels typical of the active season. Inferences of upstream regulatory molecules potentially driving these postfeeding responses implicate peroxisome proliferator-activated receptor gamma (PPARG) and other known regulators of insulin sensitivity, providing new insight into high-level regulatory mechanisms involved in shifting metabolic phenotypes between hibernation and active states.

Fitness consequences and ancestry loss in the Apennine brown bear after a simulated genetic rescue intervention.

Reduction in population size, with its predicted effects on population fitness, is the most alarming anthropogenic impact on endangered species. By introducing compatible individuals, genetic rescue (GR) is a promising but debated approach for reducing the genetic load unmasked by inbreeding and for restoring the fitness of declining populations. Although GR can improve genetic diversity and fitness, it can also produce loss of ancestry, hampering local adaptation, or replace with introduced variants the unique genetic pools evolved in endemic groups. We used forward genetic simulations based on empirical genomic data to assess fitness benefits and loss of ancestry risks of GR in the Apennine brown bear (Ursus arctos marsicanus). There are approximately 50 individuals of this isolated subspecies, and they have lower genetic diversity and higher inbreeding than other European brown bears, and GR has been suggested to reduce extinction risks. We compared 10 GR scenarios in which the number and genetic characteristics of migrants varied with a non-GR scenario of simple demographic increase due to nongenetic factors. The introduction of 5 individuals of higher fitness or lower levels of deleterious mutations than the target Apennine brown bear from a larger European brown bear population produced a rapid 10-20% increase in fitness in the subspecies and up to 22.4% loss of ancestry over 30 generations. Without a contemporary demographic increase, fitness started to decline again after a few generations. Doubling the population size without GR gradually increased fitness to a comparable level, but without losing ancestry, thus resulting in the best strategy for the Apennine brown bear conservation. Our results highlight the importance for management of endangered species of realistic forward simulations grounded in empirical whole-genome data.

Consecuencias en la aptitud y pérdida de ascendencia del oso pardo de los Apeninos después de un rescate genético simulado Resumen La reducción del tamaño poblacional, con los previsibles efectos sobre su aptitud, es el impacto antropogénico más alarmante sobre las especies amenazadas. Mediante la introducción de individuos compatibles, el rescate genético (RG) es una estrategia prometedora para reducir la carga genética revelada por la endogamia y restaurar la aptitud de las poblaciones en declive, aunque todavía se debate la eficiencia de esta. Aunque el RG puede mejorar la diversidad genética y la aptitud, también puede producir pérdida de ascendencia, lo que puede dificultar la adaptación local, o sustituir con variantes introducidas por los migrantes los acervos genéticos únicos que han evolucionado en grupos endémicos. En este trabajo realizamos simulaciones genéticas a futuro basadas en datos genómicos empíricos para evaluar los beneficios del RG en términos de aptitud y los riesgos de la pérdida de ascendencia en el oso pardo de los Apeninos (Ursus arctos marsicanus). Quedan aproximadamente 50 individuos de esta subespecie aislada que cuentan con una menor diversidad genética y un mayor nivel de endogamia comparado con otros osos pardos europeos y se ha sugerido que el RG podria reducir el riesgo de extinción de esta población. Comparamos 10 escenarios de RG en los que variaban el número y las características genéticas de los osos migrantes con un escenario sin RG con aumento demográfico causado por factores no genéticos. La introducción de 5 individuos procedentes de una población europea de oso pardo con mayor aptitud o niveles menores de mutaciones deletéreas que el oso pardo de los Apeninos produjo un rápido aumento de la aptitud del 10-20% en la subespecie y hasta un 22.4% de pérdida de ascendencia durante 30 generaciones. En las simulaciones sin un aumento demográfico, la aptitud empezó a disminuir de nuevo después de unas pocas generaciones. La duplicación del tamaño de la población sin RG aumentó gradualmente la aptitud hasta un nivel comparable al de algunos escenarios de RG, pero sin pérdida de ascendencia, por lo que parece ser la mejor estrategia para la conservación del oso pardo de los Apeninos. Nuestros resultados resaltan la importancia que tienen las simulaciones realistas a futuro basadas en datos empíricos del genoma completo para la gestión de especies amenazadas.

Age-related alterations in metabolome and microbiome provide insights in dietary transition in giant pandas.

We conducted UPLC-MS-based metabolomics, 16S rRNA, and metagenome sequencing on the fecal samples of 44 captive giant pandas (Ailuropoda melanoleuca) from four age groups (i.e., Cub, Young, Adult, and Old) to comprehensively understand age-related changes in the metabolism and gut microbiota of giant pandas. We characterized the metabolite profiles of giant pandas based on 1,376 identified metabolites, with 152 significantly differential metabolites (SDMs) found across the age groups. We found that the metabolites and the composition/function of the gut microbiota changed in response to the transition from a milk-dominant diet in panda cubs to a bamboo-specific diet in young and adult pandas. Lipid metabolites such as choline and hippuric acid were enriched in the Cub group, and many plant secondary metabolites were significantly higher in the Young and Adult groups, while oxidative stress and inflammatory related metabolites were only found in the Old group. However, there was a decrease in the α-diversity of gut microbiota in adult and old pandas, who exclusively consume bamboo. The abundance of bacteria related to the digestion of cellulose-rich food, such as Firmicutes, Streptococcus, and Clostridium, significantly increased from the Cub to the Adult group, while the abundance of beneficial bacteria such as Faecalibacterium, Sarcina, and Blautia significantly decreased. Notably, several potential pathogenic bacteria had relatively high abundances, especially in the Young group. Metagenomic analysis identified 277 CAZyme genes including cellulose degrading genes, and seven of the CAZymes had abundances that significantly differed between age groups. We also identified 237 antibiotic resistance genes (ARGs) whose number and diversity increased with age. We also found a significant positive correlation between the abundance of bile acids and gut bacteria, especially Lactobacillus and Bifidobacterium. Our results from metabolome, 16S rRNA, and metagenome data highlight the important role of the gut microbiota-bile acid axis in the regulation of age-related metabolism and provide new insights into the lipid metabolism of giant pandas. IMPORTANCE The giant panda is a member of the order Carnivora but is entirely herbivorous. The giant panda's specialized diet and related metabolic mechanisms have not been fully understood. It is therefore crucial to investigate the dynamic changes in metabolites as giant pandas grow and physiologically adapt to their herbivorous diet. This study conducted UPLC-MS-based metabolomics 16S rRNA, and metagenome sequencing on the fecal samples of captive giant pandas from four age groups. We found that metabolites and the composition/function of gut microbiota changed in response to the transition from a milk-dominant diet in cubs to a bamboo-specific diet in young and adult pandas. The metabolome, 16S rRNA, and metagenome results highlight that the gut microbiota-bile acid axis has an important role in the regulation of age-related metabolism, and our study provides new insights into the lipid metabolism of giant pandas.

Formerly bile-farmed bears as a model of accelerated ageing.

Bear bile-farming is common in East and Southeast Asia and this farming practice often results in irreversible health outcomes for the animals. We studied long-term effects of chronic bacterial and sterile hepatobiliary inflammation in 42 Asiatic black bears (Ursus thibetanus) rescued from Vietnamese bile farms. The bears were examined under anesthesia at least twice as part of essential medical interventions. All bears were diagnosed with chronic low-grade sterile or bacterial hepatobiliary inflammation along with pathologies from other systems. Our main finding was that the chronic low-grade inflammatory environment associated with bile extraction in conjunction with the suboptimal living conditions on the farms promoted and accelerated the development of age-related pathologies such as chronic kidney disease, obese sarcopenia, cardiovascular remodeling, and degenerative joint disease. Through a biomimetic approach, we identified similarities with inflammation related to premature aging in humans and found significant deviations from the healthy ursid phenotype. The pathological parallels with inflammageing and immuno-senescence induced conditions in humans suggest that bile-farmed bears may serve as animal models to investigate pathophysiology and deleterious effects of lifestyle-related diseases.

A multi-tissue gene expression dataset for hibernating brown bears.

OBJECTIVES: Complex physiological adaptations often involve the coordination of molecular responses across multiple tissues. Establishing transcriptomic resources for non-traditional model organisms with phenotypes of interest can provide a foundation for understanding the genomic basis of these phenotypes, and the degree to which these resemble, or contrast, those of traditional model organisms. Here, we present a one-of-a-kind gene expression dataset generated from multiple tissues of two hibernating brown bears (Ursus arctos). DATA DESCRIPTION: This dataset is comprised of 26 samples collected from 13 tissues of two hibernating brown bears. These samples were collected opportunistically and are typically not possible to attain, resulting in a highly unique and valuable gene expression dataset. In combination with previously published datasets, this new transcriptomic resource will facilitate detailed investigation of hibernation physiology in bears, and the potential to translate aspects of this biology to treat human disease.